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Procede et dispositif de detection de la fibrillation auriculaire cardiaque par 
la methode des ondelettes. 

DOMAINE TECHNIQUE ET ART ANTERffiUR 
5 L'invention conceme le domaine de Vanalyse et de Interpretation des 

variations de iythme se produisant dans un signal ou un phenomene physique qui 
est initialement regulier, puis qui devient irregulier, et qui peut ensuite redevenir 
regulier. 

Elle s'applique en particulier a I'analyse des variations du rythme cardiaque. 
10 Elie concerne en particulier un procede et un dispositif pour d6tecter un 

trouble du rythme cardiaque appele fibrillation auriculaire (FA). 

La fibrillation auriculaire (FA) cbnstitue le trouble du rythme le plus 
frequemment rencontre en pratique* cardiplpgique; . ■-:•!.• 

Dans un electrocardiogramme, et comme illustre sur la figure 1, on appelle 
15 intervalle RR l'intervalle de temps separant deux maxima 2, 4 de complexes QRS 
d*un electrocardiogramme. 

Les arythmies cardiaques supraventriculaires (FA) peuvent amener a de 
s^veres probtemes fonctionnels et etre associees a une grande mortalite 
cardiovasculaire. Le nombred'arythmieaugmente avec l'Sge. Le Cardiovascular 
20 Health Study Collaboration Research Group a montre, dans un article de B.M. 
Psaty et al. intitule "Incidence of and risks factor for atrial fibrillation in older 
adults" paru dans Circulation^ 1997, vol; 96, pages 2455-2461, au moins une 
anomalie majeure d ! electrocardiogramme (ECG) dans 29% d'un groupe de 5150 
sujets ages de 65 ans ou plus. Dans ce groupe, Tincidence de la fibrillation 
25 auriculaire (FA) augmentait, entre les groupes ages de 65-74 et 75-84 ahs, de 1,76 
pour 100 personnes par an a 4,27, respectivement, pour les hommes, et de 1,01 a 
2, 1 6, respectivement, pour les femmes. 

D'autres auteurs (S. Rajala et al., "Prevalence of ECG findings in very old 
people", Eur Heart J, 1984, 5, 168-174) ont decrit une prevalence de la FA de 
30 l'ordre de 83% dans une population agee de 85. ans ou plus et ont aussi montre que 
le risque clinique associe augmentait avec Vage. 

Aux Etats-Unis, un accident vasculaire cerebral ischemique sur quatre, ce 
qui signifie au total 75 000 accidents par an, est associe a une FA, comme le 
montre Particle de AJ. Camm et al. "Epidemiology and mechanism of atrial 
35 fibrillation and atrial flutter", am X Cardiol. 1996, 78, Pages 3-1 1. 
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La pr&ence de FA est un puissant predicteur de I'incidence et de la 
recidive des accidents, ainsi que de Faugmentation de la mortalite apres un 
accident, a la fois sur le court et le long terme. Par rapport aux groupes etudies, le 
risque absolu d'accident attribue a la FA variait de 3% a 67% avec un risque relatif 
5 s'etalant de 1,0 a 6,9% (K.H. Fleget et al., Risk of stroke in non- rheumatic atrial 
fibrillation, Lancet, 1987, 1, pages 526-529). Le risque absolu tfinfarctus 
cerebraux cliniquement silencieux etait dans ies m&nes proportions, tandis que le 
risque relatif etait estime k 2,5 (voir a ce sujet : H. Yamanouchi et al. "Embolic 
brain infarction in nonrheumatic atrial fibrillation : a clinico-pathologic study in 

10 the elderly, Neurology, 1997, 48, pages 1593-1597). 

On pense actuellement que le type de FA, permanente ou paroxystique, ne 
modifie pas le risque d'apparition d'accident 

Toutes ces donnees font du diagnostic de la FA une priorite de sante 
majeure (D.P. Zipes : "An overview of arhythmies and antiarrhythmic 

1 5 approaches", J. Cardiovasc Electro Physiol, 1 999, 1 0, pages 267-27 1). 

Piusieurs approches ont ete considerees. La detection de la FA a ete 
realisee en utilisant des enregistrements intracardiaques au niveau auriculaire, 
aussi bien qu'au niveau du systeme de conduction ; cependant, ces techniques 
ciblaient le declenchement des defibrillateurs implantables. D f autres essais ont ete 

20 tentes sur Tactivite auriculaire moyenne. Les signaux oesophagiens pourraient 
donner le meilleur pouvoir de detection, mais sont d'une utilisation difficile sur 
une base standard. Les enregistrements ambulatoires d'ECG ont montre etre les 
plus adaptes pour detecter de tels troubles du rythme, en terme de facilite et de 
leur capacite a detecter des formes paroxystiques frequentes. 

25 Avec YECG de surface, du type de ceux obtenus par enregistrements 

ambulatoires ("Holter"), la situation idSale serait de pouvoir reconnaitre les ondes 
P (identifiees sur la figure 1), puisqu f elles sont la cle des diagnostics et de la 
classification des arythmies auriculaires ; cependant, le faible signal obtenu dans 
la plupart des cas rend cette option difficile. Dans ce contexte, piusieurs tentatives 

30 ont ete realisees pour d^velopper des methodes efficaces de reconnaissance de la 
FA dans les enregistrements Holter. Le simple graphique de Lorentz permet de 
reconnaitre des formes specifiques de la FA. D'autres methodes fiirent appliquees 
pour ameliorer la lecture de TECG, comme celle de Murgatroyd et al (F.D. 
Murgatroyd FD et al. : "Identification of atrial fibrillation episodes in ambulatory 

35 electrocardiographic recordings : validation of a method for obtaining labeled R-R 
interval files", Pace, 1995; 18 : 1315-1320). 
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Des reseaux de neurones ont 6te aussi appliques k FECG, analysant les 
intervalles RR, ainsi que les fluctuations de la ligne de base, comme decrit dans 
les articles de J.F. Yang et al. : "Artificial neural network for the diagnosis of 
atrial fibrillation", med & biol eng & comput 1994 ; 32 : 615-619 ; "Deterministic 
5 logic versus software-based artificial neural networks in the diagnosis of atrial 
fibrillation", /. Electrocardiol. 1993 ; 26 (suppl) : 90-94 ; ou dans Tarticle de 
Cubanski D. et al : "A neural network system for detection of atrial fibrillation in 
ambulatory electrocardiograms", J. Cardiovasc. Electrophysiol 1994 ; 5 : 602- 
608. 

10 Toutefois, dans les systemes Holter, ou les enregistreurs d f 6venements 

utilises communement sur des durees de Tordre de la journee, Identification de la 
FA doit toujours etre amelioree en augmentant a la fois la sensibility et la 
specificite des algorithmes de detection. En effet, les sensibilites mentionnees 
dans les articles de T.F. Yang et al. et de D. Cubanski et al. cites ci-dessus s'etaient 

15 de 82,4% a 96,6% et les specificity de 92 a 92,3%. 

EXPOSE DE LTNVENTION ■ 

L'invention a pour objet un procede de traitement d*un signal d'intervailes 
RR d'un electrocardiogramme, ce signal 6tant pr6enregistr6 sur un support de 
20 stockage de donnees ou dans des moyens de memorisation, comportant: 

- detection des zones de forte variabilite de frequence cardiaque dans le 
signal d ! intervalles RR, 

- classification de ces zones en zones de rythme sinusal (RS) et de 
fibrillation auricuiaire (FA). 

25 La detection des zones de forte variabilite de frequence cardiaque peut etre 

realisee par projection ou transformation du signal de RR dans une base ou une 
famille d'ondelettes discretes, et calcui des coefficients d(k) donnant la projection 
du signal de RR sur chaque fonction ondelette discrete. 

Par exemple, on realise une transformation en ondelettes dyadiques. 

30 Afin de selectionner les coefficients superieurs a ime valeur TniTiimiirp, on 

procede a une etape de seuillage des coefficients dj(k) pour plusieurs niveaux j de 
decomposition. 

En outre, une etape de filtrage des coefficients seuilles peut etre prevue 
pour obtenir des zones homogenes de forte energie. 
35 Eventuellement, les variations du signal RR permettent de realiser un 

recaiage du debut et de la fin des episodes FA. 
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Les zones de forte variability de frequence cardiaque peuvent etre classees 
en zones RS et FA par calcui de Fexposant fractal de chacune de ces zones et 
comparaison dune valeur seuil Ho. 

De preference, 0,5 < H 0 < 0,9 ou encore 0,7 < Ho < 0,85. Cest dans ces 
5 plages que la sensibility et la specificity sont les meilleures. 

L'invention concerne 6galeraent un dispositif pour analyser des signaux 
numeriques dHntervaUes RR d'yiectrocardiogrammes, comportant: 

- des moyens pour, ou specialement programmes pour, dStecter des 
zones de forte variability de frequence cardiaque, 
10 - des moyens pour, ou specialement programmes pour, classifier ces 

zones en zones de rythme sinusal et en zones de fibrillation auricullaire. 

Les moyens pour detecter des zones de forte variability de frequence 
cardiaque comportent des moyens pour projeter ou transformer le signal de RR 
dans une base ou une famiile d'ondelettes discretes, et calculer des coefficients 
15 <L(k) dormant la projection du signal de RR sur chaque fonction ondelette discrete. 

Selon un autre aspect, les moyens pour classifier les zones de forte 
variability en zones de rythme sinusal et en zones de fibrillation auricullaire 
comportent des moyens pour calculer l'exposant fractal de chacune de ces zones et 
pour le comparer a une valeur seuil Ho 
20 L'invention concerne egalement un systyme de dytermination du 

rythme cardiaque, comportant un ou plusieurs enregistreurs pour enregistrer des 
electrocardiogrammes de patients, et un dispositif selon Tinvention et tei que 
decrit ci-desssus. 



25 

BREVE DESCRIPTION DES FIGURES 

Les caractyristiques et avantages de l'invention apparaitront mieux a la 
lumiere de la description qui va suivre. Cette description porte sur les exemples de 
realisation, donnes a titre explicatif et non limitatif, en se refyrant k des dessins 
30 annexes sur lesquels: 

- la figure 1 represente une partie d'un yiectrocardiogramme, avec deux 
signaux QRS, 

- la figure 2 reprysente un enregistrement Hoiter, 

- les figures 3 et 4 represented des etapes d*un procyde selon Tinvention, 
35 - la figure 5 represente des etapes pom- obtenir un signal RR k partir d f un 

signal d'ECG, 
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- la figure 6 represente des etapes de detection des zones de forte 
variability d T un signal RR, 

- la figure 7 repesente des niveaux de decomposition en ondelettes 
d'intervales RR d f un episode de fibrillation auriculaire, 

5 - la figure 8 represente la densite spectrale de puissance d*un signal de type 

RS, 

- la figure 9 represente la densite spectrale de puissance d f un signal de type 

FA, 

- la figure 10 represente un patient muni d*un systeme de mesure de type 
10 "Holter", 

- la figure 1 1 represente un enregistreur de type "Holter", 

- les figures 12 et 13 represented des moyens pour, analyser des signaux 
RR selon Tinvention, 

15 

DESCRIPTION DETAJLLEE DE MODES DE REALISATION DE 
LTNVENTION 

L'invention va d'abord etre expliquee en liaison avec les organigrammes 
20 des figures 3 a 6. 

Au prealable, on suppose qu'un signal RR a ete obtenu (etapes 3-1 et 4-1 
des figures 3 et 4), par exemple a partir d ! un electrocardiogramme (ECG) fourni a 
I'aide d'un dispositif Hotter standard. 

Un enregistrement RR ainsi obtenu est donne en figure 2. Sur cette figure, 
25 sont represents plus de 70000 intervalles RR consecutife obtenus k partir d'un 
ECG realise sur 24 heures. En ordonnee est indiquSe 1'amplitude du RR et en 
abscisse 1'indice de la serie. 

^invention s'applique aussi au traitement d'ECG obtenus sur de plus 
courtes periodes. 

30 Dans une premiere 6tape 3-2 (figure 3), on identifie les periodes d'interet, 

par detection des zones de forte variability de la frequence cardiaque. En effet, la 
fibrillation auriculaire (FA) augmente cette variability 

Une localisation precise des evenements de haute variability de frequence 
cardiaque peut etre realisee en utilisant prealablement une transformee en 

35 ondelettes discretes (TOD) du signal RR (figure 4, etape 4-2). 
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La technique de transformation en ondelettes continues est decrite dans 
Particle de S. Mallat "A theory for muitiresolution signal decomposition : the 
wavelet representation 11 , IEEE transactions on pattern analysis and machine 
intelligence, 1989, Volume 11, pages 674-693. 
5 La transformee des ondelettes discretes (TOD) est decrite dans le meme 

article de S. Mallat n s r agit d'une analyse par une familie de fonctions creee 
par les dilatations et translations d ! une ondelette mfere V|/, avec 
Vjk(x) = 2 j/2 vj/(2 j x-k), et j,k € Z (ensemble des nombres entiers relatifs). La 
projection ou la transformation de chaque fonction f (signal RR) dans la base des 

10 \|/jk est notee : Wf(j\k)= <f, \\f s ^> 9 oil <.,,> represente le produit scalaire dans 
L 2 (R). Pour un signal compose de 2 J valeurs, cette transformation am£ne k une 
serie de 2 j " 1 coefficients notes dj(k), donnant la contribution a la projection du 
signal a la position 2*k a 1'echeile 2*. La mise au carre des coefficients obtenus 
fournit la puissance du signal transforme pour chaque indice de temps et 

15 d'teheile. Differentes zones temporelles sont caracterisees par une puissance 
calculee sur les plus hautes 6chelles de leur TOD (faible j). 

La tranform^e discrete en ondelettes d'un signal discret est calculee par 
un algorithme k base de filtres miroirs en quadrature. De teis filtres sont donnes 
dans Touvrage de A.Cohen, Ondelettes et traitement numerique du signal, 

20 Edition Masson, Paris, 1 992. 

A partir du signal RR numerise, une premiere decomposition est realisee, 
qui donne un premier niveau de detail (par filtrage passe-haut) puis un premier 
niveau d'approximation (par filtrage passe-bas). 

Une deuxieme decomposition, par filtrage passe-haut et passe-bas du 

25 premier niveau d'approximation, fournit respectivement un deuxieme niveau de 
detail puis un deuxieme niveau ^approximation, qui est a son tour traite par une 
troisieme decomposition, ..etc. 

Les filtres miroirs en quadrature permettent de rdaliser les differents 
filtrages passe-haut et passe-bas successife. 

30 La figure 7 represente plusieurs niveaux de la decomposition en 

ondelettes des intervalles RR d'un episode de FA. Les abscisses correspondent a 
Tindice k et les ordonnees a Tindice j. Autrement dit, chaque ligne j donne 
Pensemble des coefficients dj(k). Au centre de cette figure une haute variability 
due a un episode de FA est identifiable. 

35 Les periodes de forte variability ayant ete identifiees, elles sont separees 

en periodes de RS et de FA (etapes 3-3 et 4-3 des figures 3 et 4). En effet, une 
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grande variability d f amplitude similaire peut aussi etre trouvee chez les sujets 
normaux, par exemple la nuit, chez certains jeunes adultes. 

A cette fin, les periodes de grande variability sont separees par analyse 
fractals 

5 Le calcul, par transferase de Fourier, de la densite spectrale de puissance 

(DSP), montre qu*un rythme sinusal normal (RS) possede une tendance generate 
en 1/f^, done presente une forme iineaire de pente negative quand elle est 
affichee en log-log, comme iliustre sur la figure 8. Sur cette figure, le pic 10 
correspond a une frequence de ventilation. La DSP de la FA, par contre, presente 

10 deux pentes differentes comme on le voit sur la figure 9 : une dans une bande de 
frequences elevees (10" 1 - 10" 2 Hz) et la seconde dans une bande de frequences 
plus basses (10°- lO^Hz). 

La pente dans la bande de haute frequence pour une FA est plus 
importante (P « 0) que la pente obtenue dans les RS (P « 1). 

15 Pour mod61iser la loi de puissance en l/f^, on applique le concept de 

mouvement Brownien fractionnaire (mBf), processus auto-similaire et non- 
stationnaire. Ce concept est decrit, par exemple, dans B. Mandelbrot et aL " 
Fractional brownian motions, fractional noises and application", SIAM Review, 
1968, volume 10, pages 422-437. 

20 La propriete d'auto-similarite est caracterisee par Pexposant de Hurst H 

qui est directement relie 4 p par la relation : p = 2H + 1. Le calcul de H est donne 
par exemple dans Particle de M Akay et aL "Fractal analysis of HRV signals : a 
comparative study " paru dans Methods of information in Medecine, 1997, 
volume 36, pages 271-273. 

25 Uincrement du processus du mouvement Brownien fractionnaire peut 

aussi etre utilise pour estimer Texposant H de Hurst, via la TOD, par : 
var [d 2j (n)] = 2 (2H ' l) «* l > a 2 (2-2 2IW ) 

ou les d 2 j sont les coefficients des details & Techelle j. On obtient alors 
liquation suivante : 
30 log 2 [var [d 2j (n)]] = (2H - l)j + f(H,a) 

ou a est Tecart-type du signal, et f une fonction non connue. 
La partie hautes-frequences d'un signal de fibrillation auriculaire 
correspond a im bruit blanc, pour lequel H « 0,5. Au contraire, un signal normal 
s f eloigne d f un bruit blanc, le coefficient de H etant superieur a 0,5. 
35 On fixe done un seuil Ho, compris entre 0,5 et 1, par exemple entre 0,5 et 

0,9 ou entre 0,7 et 0,85, par exemple H 0 = 0,8. 
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De tres bonnes valeurs de sensibilite et de sp6cificit6 de la m&hode ont 
ete obtenues pour 0,7< Ho < 0,85. En particulier, on a obtenu une sensibilite de 
99,2% pour une FA permanente et de 95,82% pour une FA paroxystique, et une 
sp6cificite de 99,96% pour un RS normal et 91,55% pour un FA paroxystique. 
5 La figure 5 represente des etapes preliminaires pour obtenir un signal RR 

a partir d'un signal d'electrocardiogramme numeris6 (ECG). Les deux premieres 
etapes (5-1 et 5-2) sont rSalisees de maniere classique dans un dispositif de type 
"Holter" standard 

Le typage permet de classifier chaque complexe en normal, ventriculaire, 
1 0 ou supraventricular. 

Les complexes QRS nonnaux et supraventriculars presentent des ondes 
P, devant le pic maximum, tandis que les anormaux (ventricuiaires) n'en 
presentent pas. Les QRS correspondant a ces derniers sont en fait des battements 
nes spontanement dans les ventricules du cceur. As sont en general plus larges 
1 5 que les nonnaux. On peut, pour les detecter, en determiner la largeur. 

Ainsi est cree un tableau ou une base de donn6es dormant le type de QRS 
pour chaque QRS, repere par son indice dans la suite des QRS enregistr6s. 

Les complexes normaux et supraventriculaires sont selectionnes a partir 
de ce tableau (etape 5-3). Un signal RR peut ensuite etre obtenu (etape 5-4) en 
20 caiculant i'ecart entre les pics des QRS successife (voir figure 1). 

Une transformee en ondelettes discretes particulierement bien adaptee 
dans le cadre de la presente invention est la transformee en ondelettes dyadiques. 

En effet, il se pose de maniere generale un probleme de reversibilite de la 
transformee en ondelettes. Si le debut d f un signal RR est correctement identifie, 
25 sa transformee en ondelettes discretes est reversible : on peut, k partir de chaque 
transformee, retrouver le signal initial par une operation d'inversion. 

Si, au contraire, I'origine du signal RR n'est pas correctement 
selectionnee, Top6ration de transformee inverse amene une distorsion et ne 
permet pas de retrouver le signal RR initial. 
30 H en resulte alors un manque de precision dans Tidentification de Torigine 

ou du debut de la fibrillation auriculaire, ce qui peut etre tres gSnant, en 
particulier pour les evenements FA de courtes durees. 

L*utilisation de la transformee dyadique, invariante par translation, 
permet de resoudre le probleme. 
35 Des informations sur ces transformees dyadiques sont donnees dans 

l'ouvrage de A. Antoniadis et G.Oppenheim, Springer Verlag, 1995, et 



WO 01/80042 PCT/FR01/01170 



notamment dans les articles suivants contenus dans cet ouvrage: "The 
Stationnary Wavelet Tranform And Some Satistical Applications" G.P.Nason et 
al., page 281-299 et "Translation Invariant De-Noising", R.R. Coifinan et aL, 
pages 125-150. 

5 La figure 6 represente des etapes pour realiser la detection ou 

Tidentification des zones de forte variability du signal RR (etape 4-3). 

Les coefficients obtenus par analyse par ondelettes sont compares aun 
seuil predetermine, afin de ne retenir que les coefficients de plus haute energie. 
Ainsi, ce seuillage peut Stre realise sur les niveaux de decomposition j = 2, 4, et 
10 8 (etape 6-1). 

Puis (etape 6-2) il est procede a un fiitrage de ces coefficients, afin 
d'obtenir des zones homogenes de forte energie. 
Pour ce faire, on procede a: 

- une recherche de debut de chaque episode ou de chaque zone de forte 
15 energie,, 

- une comptabilisation des coefficients dans Tepisode ou dans la zone, 

- une determination de fin de Tepisode ou de la zone. 

Au cours de I'etape de recherche de debut d^pisode, on compte le nombre 
de coefficients qui, dans une certaine fenetre, sont sup6rieurs au seuil 
20 predetermine. Lorsque ce nombre est superieur a un certain nombre Ns 
predetermine, il y a debut d f un episode. 

Dans la phase de comptabilisation, on compte le nombre d'evenements 
successifs qui sont inferieurs au seuil. Si ce nombre est superieur a un nombre 
predetermine, une fin d'episode est detectee. 
25 Cette etape 6-2 permet une determination assez precise des evenements. 

Dans certains cas, ce n'est cependant pas suffisant et un degre superieur 
de precision est necessaire. 

Pour cette raison, on procede a un recalage du debut et de la fin des 
Episodes (etape 6-3). 

30 Cette etape recalage permet de ne pas manquer ou laisser passer des 

evenements uniques et, de plus, d'eviter de mal identifier le debut d'un episode 

en cas d'episode plus long. 

Des evenements uniques peuvent survenir, qui sont tout aussi importants 

que les evenements successifs. Pour cette raison, ils peut etre tout aussi 
35 important, pour le praticien ou le medecin cardiologue, de disposer de 

rinformation sur les evenements simples aussi bien que sur les episodes 
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comptant plusieurs evenements. En outre, un mauvaise positionnement du debut 
de Episode, qui peut notamment apparaitre en representation sur un ecran de 
visualisation, peut troubler ce praticien ou ce medecin et le faire douter de la 
fiabilitS de Tappareil mettant en ceuvre la methode d'analyse seion invention. 
5 Pour resoudre ce probleme, apres avoir identifie les zones homogenes de 

forte energie du signal RR tranforme par ondelettes (etape 6-2), on determine ies 
zones de forte variability du signal numerique RR (6tape 6-3). 

Selon un exemple, on reprend le signal numerique RR, et on detecte les 
variations de ce signal qui sont superieures a un certain pourcentage 
10 predetermine, par exemple 20 % (on predetermine un seuil compris entre 10% et 
40% ou entre 15% et 25%). 

Lorsqu ! il y a une telle variation d ! un point a un autre du signal RR, on 
peut precisement dire qu'il y a debut d'un evenement en ce point. 

On determine de la m6me maniere la fin des episodes ou des Svenements. 
15 Ainsi, pour chaque episode, on peut obtenir un intervalle Nl - N2, qui identifie 
le rang des signaux RR compris entre debut et fin des episodes. 

Cet intervalle est ensuite reporte sur les difiterents niveaux de 
decomposition (figure 7) et seuls les coefficients dans cet intervalle sont retenus 
comme correspondant a i'episode. 
20 II est possible, ainsi, que des coefficients, consideres apres Tetape 6-2 

comme faisant partie du debut ou de la fin d'une zone homogene de forte 
energie, soient exclus de cette zone du fait d f une legere modification, due au 
recalage, du debut ou de la fin d'un episode, ou bien que des coefficients situes a 
la transition des zones homogenes de forte energie et des zones de faible energie 
25 soient reclass6s dhine zone k Tautre. 

La figure 10 represente un patient 20 muni d f un dispositif 22 de type 
"Hotter", avec ses electrodes 24,26,28,30,32 disposees sur le corps du patient. Le 
Holter enregistre les donnees d'ECG pendant ime certaine periode, par exemple 

30 environ 24 heures ou sur une periode plus courte (par exemple quelques 
minutes), ou sur des periodes plus longues (par exemple de plusieurs jours ou de 
plusieurs semaines), , et detecte et ciassifie les complexes Qrs. 

Comme illustr£ sur la figure 11, les donnees sont par exemple 
enregistrees sur une carte memoire 34 de type memoire Flash, par exemple 

35 d f environ 20 Mbytes. D'autres appareils permettent de faire des enregistrements 
sur des cassettes inserees dans le Holter 22. 
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Les figures 9 et 10 represented sch&natiquement les composants 
electroniques essentiels d'un appareil Holter. 

La figure 12 represente un dispositif 40 pour exploiter des signaux 
enregistres sur un support memoire 34, par exemple une carte m&noire. Le 
5 support 34 est introduit dans un lecteur 36, par exemple un lecteur PCMCIA, et 
les donnees subissent ensuite un traitement selon les etapes d6crites ci-dessus. 

Le dispositif 40 est done programme pour rSaliser ces Stapes de 
traitement 

Un tel dispositif comporte un micro-ordinateur PC 42 auquel les donn6es 
10 lues par le lectreur 36 sont transmises via une liaison 41 . 

Plus precisement, le micro-ordinateur 42 comporte (figure 13) un 
microprocesseur 52, un ensemble de memoires RAM 54 (pour le stockage de 
donnees), une memoire ROM 56 (pour le stockage ^instructions de programme). 
Eventuellement, une carte d'acquisition de donnees 60 met les donnees 
15 fouraies par le lecteur au format requis par le micro-ordinateur. Ces divers 
elements sont relics a un bus 62. 

Des dispositifs peripheriques (ecran ou dispositif de visualisation 64, 
souris 66, clavier 68) pennettent un dialogue interactif avec un utilisateur. En 
particulier, les moyens de visualisation (ecran) 64 pennettent de fournir & un 
20 utilisateur une indication visuelle relative a PECG enregistre, aux intervalles RR 
calcules, aux coefficients calcules par transforlmee par ondelettes (comme sur la 
figure 7). Par selection dans des menus et en cliquant sur des boutons ou des 
icones, un operateur, par exemple un medecin cardiologue, peut visualiser ces 
dififerentes courbes, en entier ou par portion, avec agrandissement. Sur ces courbes 
25 peuvent apparaitre les intervalles determines comme etant des episodes de 
fibrillation auriculaire, et notamment le debut et la fin de chaque Episode. 

Dans le micro-ordinateur 42, sont chargees les donnees ou les instructions 
pour mettre en ceuvre un traitement des donnees selon Tinvention. 

Ces donnees ou instructions pour le traitement des donnees peuvent Stre 
30 transferees dans une zone memoire du micro-ordinateur 42 & partir d f une disquette 
ou de tout autre support pouvant etre lu par un micro-ordinateur ou un ordinateur 
(par exemple: disque dur, memoire morte ROM, memoire vive dynamique DRAM 
ou tout autre type de memoire RAM, disque optique compact, 616ment de stockage 
magnetique ou optique). 
35 L'afBchage des episodes classifies comme etant des episodes de 

fibrillation auricullaire ne constitue pas, en soi, un diagnostic. Seul un medecin 
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cardiologue peut, k la vue de r61ectrocardiogramme initial (qui a lui aussi et6 
numerise, mis en memoire sur le support 34 et transfere en memoire du 
microordinateur 42) confirmer ['existence d'une pathologic En particulier, le 
medecin doit encore exercer une observation de la ligne de base 3 de 

5 relectrocardiogramme (voir figure 1), afin de confirmer si un episode, 
apparemment classifie comme 6tant de type FA, est bien le signe d*une 
fibrillation auriculaire. 

Le precede de traitement presente ci-dessus est done une aide au 
diagnostic, mais pas le diagnostic lui-meme, qui demande une etape 

10 ^interpretation mSdicale supplemental. 

EXEMPLE 

Cet exemple concerne une population de 50 sujets : 19 soufifrant de 
15 fibrillation auriculaire chronique (FAC), 15 de fibrillation paroxystique (FAP), et 
16 patients controles avec un rythme sinusal normal (RS). 

Tous les sujets ont ete enregistres pendant 24 heures en utilisant un 
systeme Holter ECG numerique (DuoHoiter, Novacor, Paris, France). Les 
enregistrements ont ete directement numerises a une frequence de 100 Hz sur 
20 Tenregistreur ; ils ont ensuite 6t6 revus et edites manuellement, si necessaire, 
pour obtenir une liste d'intervalles RR consecutifs. 

Le traitement du signal a ete realise avec le iogiciel MATLAB® 5.0 et 
Vanalyse statistique des resultats avec Statview ® sur un Power Macintosh. ® 
Les resultats sont presentes sous la forme moyenne±6cart-type. 
25 Les resultats sont les suivants : 

a. Enregistrements ECG du groupe controle avec rythme normal 
L'examen standard des enregistrements ECG sur 24 heures des 16 sujets 
RS a montre Tabsence de toutes arythmie cardiaque. 

Chez un sujet, deux episodes de FA furent faussement detectes (duree de 
30 154 et 193 battements respectivement ; frequence cardiaque moyenne de 
81,5±14,8 et 76,0±9,5 battements par minutes (bpm) respectivement). II n'y avait 
pas de difference significative dans les indices temporels entre le sujet avec les 
fausses detections et les 15 autres. 

Ainsi, la specificite de la detection de la FA basee sur les complexes QRS 
35 atteint 99,96% (1 083 537 Qrs vrais negatifs contre 447 faux positifs). Comme 
tous les faux positifs fiirent obtenus chez un seul patient, la specificite etait de 
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93,7% si Ton considere le patient comme unite d'analyse (16 vrais negatifs contre 
1 faux positif). 

b. Groupe des FA chr oniques 

L'examen standard des enregistrements ECG de 24 hemes a montre une 
5 FA permanente chez les 19 sujets. Les enregistrements etaient d'une duree de 
18,8±4,6 heures. Chez les 19 patients, ie nombre de QRS etait de 87953±27491 
et la frequence cardiaque de 75,4+1,0 bpm (intervalle de 23,1 k 289,7 battements 
par minute). 

Un RS a ete faussement identifie chez 7 patients sur 40 episodes 
10 (frequence cardiaque 82,5±9,5, duree 374±360 battements) ; un patient 
regroupait 24 des 40 episodes alors que les autres apparaissaient chez 6 autres 
patients- Ces episodes faux negatifs etait dus a du flutter auriculaire et de la 
tachycardie auriculaire dans, respectivement, 34 et 1 des 40 cas. Chez 5 des 7 
patients, la FA n'etait pas reconnue au debut des enregistrements ; cependant le 
1 5 delai de reconnaissance 6tait limits k 1 78,4±70 battements. 

Ainsi, la sensibilite de detection de la FA atteignait 99,2% en considerant 
les complexes QRS (2064197 QRS vrais positifs contre 16018 faux negatifc). A 
i'exception des episodes non detectes, tous les patients ont ete reconnus comme 
souffrant de FA et la sensibilit6 obtenue en considerant les patients etait de 
20 100%. La plus grande periode non detectee atteignait 2075 battements et la plus 
courte 120. Le patient avec 24 episodes non detectes avait 10490 battements en 
faux negatifs ; pour cet enregistrement particulier, la sensibilite de detection de 
la FA basee sur les complexe QRS etait de 89,9%. 
c Groupe des FA intermittentes. 
25 Les 15 sujets avec FAP avaient un total de 36 Episodes de FAP dont la 

duree etait de 2,5±3,4 heures (de 34 s a 15h36). Pendant ces episodes, la frequence 
cardiaque moyenne Stait de 86,2±20,7 (de 16,4 a 163,1) bpm et le nombre de QRS 
de 15055±22125 (de 52 k 94990 QRS). 

Seul un des 36 episodes restait non detecte (260 QRS, frequence cardiaque 
30 de 123,2±28,1 bpm). Pour les episodes detects, il y avait des deiais ou des 
avances dans la reconnaissance de la FA au debut et a la fin des 6v£nements : 

- debut des Episodes : dans 8 episodes il y avait un delai de 9U95 
battements (de 2 a 255) ; 23 episodes Stait d&ectes prematurement (dus k des 
tachycardies supraventriculaires) de 186±621 battements (de 2 a 2994). 
35 - fin des episodes : deux episodes 6tait prematures (un de 3486 battements, 

Pautre de 333) ; 21 avaient un delai de 60*103 battements (de 10 k 429). 
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Durant les episodes de FAP, 14 fausses interruptions etaient observes, 
parmi lesquelles 13 chez le meme patient, avec une duree de 1484±151 
battements ; la derniere 6tait de 2001 battements chez un autre patient ; le total 
de complexes QRS faux negatifs atteignait alors 28695 (1,9%). Ces episodes 
5 faux negatifs etait dfls a du flutter auriculaire et k de la tachycardie auriculaire 
dans, respectivement, 5 et 8 des 14 cas. 

n y avait aussi 37 fausses detections de 1573±2309 battements (de 147 a 
11491). 

Ainsi dans le groupe des FAP, la sensibilite basee sur les complexes QRS 
10 atteignait 95,8% et la specificite 91,5%/ 

Tous les patients du groupe etaient reconnus comme ayant de la FA, done 
la sensibilite basee sur les patients etait de 100%. 

La transformee en ondelettes et Tanalyse fractale associees selon 
i'invention ont fourni xine sensibilite de 95% pour detecter les complexes QRS en 
15 FA, et une specificite de 99,6%, avec une limitation due aux delais de 
reconnaissance de la FA au debut et h la fin des episodes d'arythmie. En se basant 
sur les patients, la specificite atteignait 93% dans le groupe controle normal et la 
sensibilite 100% que ce soit dans le groupe FAC ou dans le groupe FAP. 

Les enregistrements intra-auriculaires ont souvent ete utilises pour 
20 detecter la FA, ainsi que ['analyse de Pactivit6 electrique de i'oreiilette au niveau 
de Foesophage, avec une sensibilite et une sp6cificit6 s'etalant de 52,4-96,6% et 
92,0-92,3%, respectivement (voir notamment S. Rajaia et al., "ECG findings and 
survival in very old people", Eur Heart J. 1985, volume 6, pages 274-252). 

De teis enregistrements n ! ont pas ete utilises etant donne que le but 6tait la 
25 detection de la FA de fa9on non-invasive par enregistrement Holter standard. 

Pour obtenir une detection automatique, des enregistrements Holter ont ete 
analyses par des reseaux de neurones (T.F. Yang et aL "Artificial neural network 
for the diagnosis of atrial fibrillation", med. & bio. Eng.& Comp. 1994, 32, pages 
615-619), dormant une sensibilite de 92% et une specificity de 92,3% dans la 
30 reconnaissance des QRS. En utilisant un reseau de neurones aliments i la fois par 
Tinformation des intervalles RR et une analyse morphologique de la ligne de base 
du trace ECG, Cubanski D. et al (J. Cardiovasc. Electrophysiol. 1994, 5, pages 
602-608), travaillant sur des groupes de 10 QRS consecutifs, ont obtenu une 
sensibilite de 82,4% et une specificite de 96,6%. Ainsi, la sensibilite et la 
35 specificite obtenues par le proced£ selon Tinvention sont meilleurs que ceiles des 
etudes precedentes. 
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Notre m&hode possede certaines limitations. I/utilisation de r analyse de 
la VFC seule est sensible a tout ce qui peut dSgrader la VFC, notamment ies 
extrasystoles supraventriculaires, les tachycardies supraventriculaires et les autres 
types d'arythmie. Les detections de faux positifs peuvent etre dues a la presence 
5 de quelques arythmies supraventriculaires, isolees ou groupees, generant une VFC 
importante ; les detections de faux positifs se sont produites sur 0,04% des QRS 
du groupe controle normal et dans 4,7% des QRS du groupe PAF. Ces QRS faux 
positifs peuvent etre reclassifies correctement en utilisant un trace enregistr6 par 
un stockage automatique d'6venements. A Pinverse, les non-detections peuvent se 
10 produire dans d'autres arythmies auriculaires, comme la tachycardie 
supraventricuiaire paroxystique, ou le flutter auriculaire, du a la relative regularity 
de ces rythmes. 

Dans des conditions cliniques, une non-d6tection est plus problematique 
qu'une detection de faux positifs, car dans ce cas le trace n'est pas stocke pour 

15 i'utilisateur, ce qui empeche toute correction ultSrieure. Dans les cas etudies, les 
non-d&ections n'6taient pas frequentes et representaient moins de 3% des 
episodes (1 sur 36) dans le groupe de FAP et moins de 2% des QRS. Des 
ameliorations de Talgorithme pourraient diminuer ces non-d&ections, au 
detriment toutefois de la specificity 

20 ^invention est done adaptee a la detection automatique de la FA dans de 

larges populations en utilisant des enregistrements Holter standard Ceci est 
particulierement important chez les personnes agees, chez lesquelles la FA, aussi 
bien paroxystique que chronique, est fr6quente et presente un risque majeur 
d'accident vasculaire cerebral. Implemente dans des appareils Holter standard ou 

25 dans des enregistreurs d'evenement a long-terme, Mgorithme ou le procede selon 
I'invention permet de mieux mettre en evidence la FA afin de prendre des mesures 
preventives avant Tapparition de complications dramatiques. 
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REVINDICATIONS 

I. Procede de traitement dhin signal d'intervalles RR d f un 6lectrocardiogramme, 
comportant: 

5 - la detection des zones de forte variability de frequence cardiaque dans le 

signal RR, 

- la classification de ces zones en zones de rythme sinusal (RS) et de 
fibrillation auriculaire (FA). 

10 2. Procede selon la revendication 1, la detection des zones de forte variability de 
frequence cardiaque etant reaiisee par projection ou transformation du signal 
d'intervalles RR dans une base ou une famille d'ondelettes discretes, et calcul des 
coefficients &(k) donnant la projection du signal de RR sur chaque fonction 
ondelette discrete. 

15 

3. Procede selon la revendication 2, la projection ou la transformation du signal 
d'intervalles RR etant reaiisee dans une base ou une famille d'ondeiettes 
dyadiques. 

20 4. Procede selon Tune des revendications 2 ou 3, comportant une etape de 
seuillage des coefficients dj(k) sur un ou piusieurs niveaux j de decomposition 
afin de sSlectionner ies coefficients superieurs a une valeur minimum. 

5. Procede selon la revendication 4, comportant en outre une etape pour obtenir 
25 des zones homogenes de forte 6nergie. 

6. Procede selon la revendication 5, l'etape pour obtenir des zones homogenes de 
forte energie comportant un filtrage des coefficients seuill6s. 

30 7. Procede selon Tune des revendications 4 & 6, comportant, pour chaque zone 
homogene de forte eneigie, une recherche du debut de la zone, une 
comptabilisation des coefficients dans cette zone, et une determination de la fin de 
la zone. 

35 8. Procede selon la revendication 7, la recherche de d6but de la zone comportant 
un comptage des coefficients qui sont superieurs au seuil predetermine. 
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9. Procede seion la revendicatioa 7 ou 8, la comptabilisation des coefficients 
comportant un comptage du nombre d'6venements successifc qui sont inferieurs a 
un seuil predetermine. 

5 

10. Proced6 selon la revendication 9, une fin d'episode etant detectee si le nombre 
d'evenements comptes est superieur a un nombre predetermine. 

11. Procede selon Tune des revendications 5 k 10, comportant en outre une etape 
10 de recalage ou ^identification des debuts et fins des episodes, a Taide des 

variations du signal RR 

12. Procede selon la revendication 11, les coefficients donnant la projection du 
signal de RR etant retenus comme faisant partie dhine zone homogene de forte 

15 energie seulement si Us sont compris entre un debut et une fin d f un episode. 

13. Procede selon la revendication 12, comportant en outre un affichage, sur des 
moyens de visualisation, d'au moins une partie de relectrocardiogramme et ime 
identification, sur cet affichage, des zones de fibrillation auriculaire. 

20 

14. ProcedS selon ftine des revendications 1 a 13, les zones de forte variability de 
frequence cardiaque etant ciassees en zones RS et en zone- FA par calcul de 
Texposant fractal de chacune de ces zones et comparaison a une vaieur seuil Ho. 

25 15. Procede selon la revendication 14, la vaieur seuil Ho etant comprise entre 0,5 
et0,9. 

16. Proced6 selon la revendication 15, la vaieur seuil 6tant comprise entre 0,7 et 
0,85. 

30 

17. Procede selon Tune des revendications 1 a 16, relectrocardiogramme ayant et6 
acquis pendant une duree comprise entre quelques minutes et quelques mois. 



18. Dispositif pour analyser des signaux num&iques d'intervalles RR 
35 d'electrocardiogrammes, comportant: 
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- des moyens pour classifier ces zones en zones de rythme sinusal (RS) et 
en zones de fibrillation auriculaire (FA). 

19. Dispositif selon la revendication 18, les moyens pour detecter des zones de 
5 forte variability de frequence cardiaque comportant des moyens pour projeter ou 
transformer le signal d ! intervalles RR dans une base ou une famille d'ondelettes 
discretes, et calculer des coefficients d.(k) dormant la projection du signal de RR 
sur chaque fonction ondelette discrete. 

10 20. Dispositif selon la revendication 18 ou 19, les moyens pour classifier les zones 
de forte variability en zones de rythme sinusal et en zones de fibrillation auriculaire 
comportant des moyens pour calculer 1'exposant fractal de chacune de ces zones et 
pour le comparer a une valeur seuil Ho. 

15 21. Dispositif selon Tune des revendications 18 a 20, comportant en outre des 
moyens de visualisation pour afficher les zones de fibrillation auriculaire. 

22. Dispositif selon Tune des revendications 18 a 20, comportant en outre des 
moyens de visualisation pour visualiser les zones de fibrillation auriculaire et les 
20 electrocardiogrammes, et pour visualiser les zones de fibrillation auriculaire par 
rapport aux electrocardiogrammes. 



25 



23. Systeme de determination du rythme cardiaque, comportant un ou plusieurs 
enregistreurs pour enregistrer des electrocardiogrammes de patients, et un 
dispositif selon Tune des revendications 18 a 22. 
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Process and device for detecting cardiac atrial fibrillation using the wavelet method. 

TECHNICAL SCOPE AND PRIOR ART 

The invention involves the areas of analysis and interpretations of variations in rhythm produced 
in a signal or physical phenomenon which is initially regular, but which then becomes irregular, and 
which may then become regular once again. 

In particular, it applies to the analysis of variations in cardiac rhythm. 

It involves in particular a process and device for detecting a cardiac rhythm problem called atrial 
fibrillation (AF). 

Atrial fibrillation (AF) represents the rhythm problem most frequently encountered in cardiology 
practice. 

In an electrocardiogram, and as illustrated in figure 1, we call the time interval separating the two 
max points 2 and 4 of QRS electrocardiogram complex functions, the RR interval. 

Supraventricular cardiac arrhythmias (AF) may lead to severe functional problems and be 
associated with a high cardiovascular mortality. The number of arrythmias increases with age. The 
Cardiovascular Health Study Collaboration Research Group has shown, in an article by B.M. Psaty et al. 
entitled "Incidence of and risks [sic] factor for atrial fibrillation in older adults" appearing in Circulation, 
1997, vol. 96, pages 2455-2461, at least one major electrocardiogram (EKG) anomaly in 29% of a group 
of 5150 subjects aged 65 years or more. In this group, incidence of atrial fibrillation (AF) increased, 
between the groups aged 65-74 and 75-84 years, from 1.76 out of 100 persons per year to 4.27, 
respectively, for men, and 1.01 to 2.16, respectively, for women. 

Other authors (S. Rajala et al., "Prevalence of ECG[sic[ findings in very old people", Eur Heart 7, 
1984, 5, 168-174) described a prevalence of FA on the order of 83% in a population aged 85 years or 
more, and also showed that clinical risk increases with age. 

In the United States, one cerebral vascular ischemic in four, which represents a total of 75,000 
incidents per year, is associated with an AF, as shown in the article by A.J. Camm et al. "Epidemiology 
and mechanism of atrial fibrillation and atrial flutter", am. 7. Cardiol. 996, 78, Pages 3-11. 
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The presence of AF is a powerful predictor of the incidence and recidivism of incidents, as well as 
of the increase in mortality after an incident, both in the short and long term. Relative to the groups 
studies, the absolute risk of incident attributed to AF varied from 3 to 67% with a relative risk ranging 
from 1.0 to 6.9% (K.H. Flegel et al., Risk of stroke in non-rheumatic atrial fibrillation, Lancet, 1987, 1, 
pages 526-529). Absolute risk of clinically silent cerebral infarctions was in the same proportions, 
whereas relative risk was estimated at 1.5 (on this subject, see: H. Yamanouchi et al. "Embolic brain 
infarction in non-rheumatic atrial fibrillation: a clinico-pathologic study in the elderly, Neurology, 1997, 
48, pages 1593-1597). 

At present we believe that the type of AF, permanent or paroxystic, does not change the risk an 
incident will occur. 

All this data make AF diagnosis a major health priority (D.P. Zipes: "An overview of arrhythmias 
and anti-arrhythmic approaches", 7. Cardiovasc Electro Physiol 1999, 10, pages 267-271). 

Several approaches were considered. Detection of AF was performed using intracardiac 
recordings at the atrial level, as well as at the level of the conduction system; still, these techniques 
targeted the tripping of implantable defibrillators. Other tests were attempted on the mean atrial activity. 
Esophageal signals would give better detection power, but are difficult to use on a standard basis. 
Ambulatory EKG recording were shown to be the most adapted for detecting this type of rhythm 
disturbances, in terms of ease and their ability to detect frequent paroxystic forms. 

With the surface EKG, of the type obtained through ambulatory recordings ("Holter"), the ideal 
situation would be to be able to recognized P waves (identified in figure 1), because they are the key to 
diagnosis and classification of atrial arrhythmias; still, the weak signal obtained in most cases makes this 
option difficult. In this context, several attempts were carried out to develop effective methods of 
recognizing AF in Holter recordings. The simple Lorentz graph allows us to recognize specific forms of 
AF. Other methods were applied to improve EKG reading, like that in Murgatroyd et al. (F.D. 
Murgatroyd FD et al.: "Identification of atrial fibrillation episodes in ambulatory electrocardiographic 
recordings : validation of a method for obtaining labeled R-R interval files", Pace, 1995; 18 : 1315-1320). 
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Neural networks were also applied to the EKG, analyzing RR intervals, as well as fluctuations of 
the base line, as described in the articles by J.F. Yang et al.: "Artificial neural network for the diagnosis of 
atrial fibrillation", med & biol eng & comput 1994; 32 : 615-619; "Deterministic logic versus software- 
based artificial neural networks in the diagnosis of atrial fibrillation", /. Electrocardiol. 1993; 26 (suppl): 
90-94; or in the article by Cubanski D. et al.: "A neural network system for detection of atrial fibrillation 
in ambulatory electrocardiograms", /. Cardiovasc. Electrophysiol 1994; 5 : 602-608. 

Even so, in Holter systems, where event recorders commonly used over the course of a day, 
identification of the AF must still be improved by increasing both sensitivity and specificity of detection 
algorithms. En effect, the sensitivities indicated in the articles by T.F. Yang et al. and D. Cubanski et al. 
Cited above range from 82.4% to 96.6% and specificities from 92 to 92.3%. 

BRIEF DESCRIPTION OF THE INVENTION 

The subject of the invention is handling of a signal, of RR interval, of an electrocardiogram, where 
this signal was pre-recorded on a data storage support or in a means for memorization, including: 

- detection of areas of high cardiac frequency variability in the RR interval signal, 

- classification of these areas into areas of sine rhythm (SR) and atrial fibrillation (AF). 
Detection of areas of high cardiac frequency variability may be carried out by projection or 

transformation of the RR signal in a base or a family of discreet wavelets, and calculation of coefficients 
d.(k) that give the RR signal projection on each discreet wavelet function. 
For example, we carried out a transformation in dyadic wavelets. 

In order to choose coefficients that are greater than a minimum value, we proceed with a threshold 
stage for coefficients dj(k) for multiple breakdown levels j. 

Moreover, a filtering stage for threshold coefficients may be provided to obtain homogeneous 
high-energy zones. 

Variations in RR signals may possibly allow us to carry out a re-setting of the start and end of AF 
episodes. 
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Zones of high variability in cardiac frequency may be classed as zones SR and AF by calculating 
the fractal exponent of each of these zones and comparing them to a threshold value Ho. 

Preferably, 0.5 < H 0 < 0.9 or even 0.7 < H 0 < 0.85. Within these ranges, sensitivity and specificity 
are the best. 

The invention also involves a device for analyzing numerical signals of electrocardiogram RR 
intervals, including: 

- a means for, or especially programmed for, detecting areas of high cardiac frequency 
variability, 

- a means for, or especially programmed for, classifying these areas into sinus rhythm areas and 
atrial fibrillation areas. 

The means for detecting areas of high cardiac frequency variability consist of a means for 
projecting or transforming the RR signal into a base or a family of discreet wavelets, and calculating 
coefficients d.(k) giving the RR signal projection on each discreet wavelet function. 

According to another aspect, the means for classifying the areas of high variability into sine 
rhythm areas and atrial fibrillation areas consist of a means for calculating the fractal exponent of each 
one of these areas and comparing it to a threshold value Ho. 

The invention also concerns a system for determining cardiac rhythm, consisting of one or more 
recorders for recording patients' electrocardiograms, and a device according to the invention and as 
described above. 

BRIEF DESCRIPTION OF FIGURES 

The characteristics and advantages of the invention will become more apparent in light of the 
description that follows. This description is supported by production examples, data that is given for 
clarification purposes and which is by no means limiting, by referring to the attached drawings in which: 

- Figure 1 represents a part of an electrocardiogram, with two QRS signals, 

- Figure 2 represents a Holter recording, 

- Figures 3 and 4 represent stages in a procedure according to the invention, 

- Figure 5 shows stages in obtaining an RR signal using an EKG signal, 
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- Figure 6 shows stages in detection of areas of high variability in an RR signal, 

- Figure 7 shows breakdown levels of RR interval wavelets in an atrial fibrillation episode, 

- Figure 8 shows the spectral power density of an SR type signal, 

- Figure 9 shows the spectral power density of an AF type signal, 

- Figure 10 shows a patient outfitted with a "Holter"-type measurement system, 

- Figure 1 1 shows a "Holter'Mype recorder, 

- Figures 12 and 13 show a means for analyzing RR signals according to the invention. 

DETAILED DESCRIPTION OF METHODS OF EMBODIMENT OF THE INVENTION 

The invention will first of all be explained in conjunction with the flow charts in figures 3 to 6. 

First of all, we suppose that an RR signal has been obtained (stages 3-1 and 4-1 in figures 3 and 
4), for example using an electrocardiogram (EKG) supplied with the help of a standard Holter device. 

An RR recording obtained in this way is given in figure 2. On this figure, we show more than 
70,000 consecutive RR intervals obtained using an EKG carried out over 24 hours. The RR amplitude is 
shown on the ordinate and the series index on the abscissa. 

The invention can also be applied to the handling of EKG's obtained over shorter periods. 

In a first stage 3-2 (figure 3), we identify the periods of interest, by detecting areas of high cardiac 
frequency variability. In effect, atrial fibrillation (AF) increases this variability. 

A precise localization of high cardiac frequency variability events may be carried out by using, 
beforehand, a transformation into discreet wavelets (TDW) of an RR signal (figure 4, stage 4-2). 
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The continuous wavelet transform technique is described in the article by S. Mallat "A theory for 
multiresolution signal breakdown : the wavelet representation", IEEE transactions on pattern analysis 
and machine intelligence, 1989, Volume 11, pages 674-693. 

Transformation of discreet wavelets (TDW) is described in thee article by S. Mallat. It is about an 
analysis by a family of functions created by dilations and translation of a mother wavelet *F, with 
¥fk(x) = 2 j/2 ¥(2 j x-k), and j, k = Z (all related whole numbers). The projection or transformation of each 
function f (RR signal) in the base is noted : W f (j,k) = <F, where <....> represents the scalar 
product in L 2 (R). For a signal composed of 2 f values, this transformation leads to a series of 2 yX 
coefficients notes df(k), giving the contribution to the projection of the signal to the position 2 J k in scale 
2K Squaring the coefficients obtained provides the power of the signal transformed for each time and 
scale index. Different temporal zones are distinguished by a power calculated over the highest scales of 
their TDW (lowj). 

Discrete transformation into wavelets of a discrete signal is calculated by an algorithm based on 
squared mirror filters. Such filters are given in the work by A. Cohen, Wavelengths and numerical 
handling of the signal, Edition Masson, Paris, 1992. 

Starting with the digitized RR signal, a primary breakdown is carried out, which gives a first level 
of detail (by low-pass filtering), than a first level of approximation (by low-pass filtering). 

A second breakdown, by high-pass and low-pass filtering of the first level of approximation, 
provides, respectively, a second level of detail, then a second level of approximation, which is, in its turn, 
handled by a third breakdown, etc. 

The squared filter mirrors allow us to carry out various successive high-pass and low-pass 
filterings. 

Figure 7 shows several levels of the breakdown into wavelets of the RR intervals of an AF 
episode. The abscissas correspond to index k and ordinates to index j. In other words, each line j gives 
all coefficients dj(k). At the center of this figure a high variability due to an AF episode can be identified. 

Where periods of high variability have been identified, they are separated into periods of SR and 
AF (stages 3-3 and 4-3 in figures 3 and 4). In effect, a 
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great variability in similar amplitude can also be found in normal subjects, for example at night, in some 
young adults. 

To this end, periods of great variability are separated by fractal analysis. 

Calculation, using Fourier transforms, of the spectral power density (SPD), shows that a normal 
sine rhythm (SR) has a general tendency in l/i*, and therefore presents a linear shape of negative slope 
when it is plotted as log-log, as shown in figure 8. In this figure, peak 10 corresponds to a ventilation 
frequency. The SPD of the AF, on the other hand, shows two different slopes, as can be seen in figure 9 : 
in a high frequency band (10 1 - 10 2 Hz) and the second in a band of lower frequencies (10" 3 - 10" 4 Hz). 

The slope in the high-frequency band for AF is greater (P « 0) than the slope obtained in the SR(P 

In order to model the law of power as l/f* 3 , we apply the concept of fractional Brownian movement 
(fBm), a process that is auto-similar and non-stationary. This concept is described, for example, in b. 
Mandelbrot et al. "Fractional Brownian motions, fractional noises and application", SIAM Review, 1968, 
volume 10, pages 422-437. 

The auto-similarity property is distinguished by the H Hurst exponent which is directly related to p 
by the relation : p - 2H + 1. Calculation of H is given, for example, in the article by M. Akay et al. 
"Fractal analysis of HRV signals: a comparative study" which appeared in Methods of information in 
Medicine, 1997, volume 36, pages 271-273. 

The increment of the fractional Brownian movement process can also be used to estimate the 
Hurst exponent H, using the TDW, by: 

Var [d 2j (n)] = 2 (2H_1)(i l) a 2 (2-2 2H1 ) 

where d2j are the coefficients of the scale details j. We thus obtain the following equation: 
log 2 [var [d 2j (n)]] = (2H - l)j + f(H,o) 

where a is the standard deviation of the signal, and f is an unknown function. 

The high-frequency part of an atrial fibrillation signal corresponds to white noise, for which H ~ 
0.5. On the other hand, a normal signal is far from white noise, since the coefficient of H is greater than 
0.5. 

Therefore we set a threshold H 0 , of between 0.5 and 1, for example between 0.5 and 0.9 or 
between 0.7 and 0.85, for example Ho = 0.8. 
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Very sensitivity and specificity values for the method were obtained for 0.7 < H 0 < 0.85. In 
particular, we obtained a sensitivity of 99.2% for a permanent AF and 95.82% for a paroxystic AF, and a 
specificity of 99.96% for a normal SR and 91.55% for a paroxystic AF. 

Figure 5 shows preliminary stages for obtaining an RR signal using a digitized electrocardiogram 
signal (EKG). The first two stages (5-1 and 5-2) are carried out in the standard manner in a device of the 
standard "Holter" type. 

Typing allows us to classify each complex function as normal, ventricular, or supraventricular. 

Normal and supraventricular QRS complex functions show P waves, before the maximum peak, 
whereas abnormal (ventricular) ones don not show them. QRS's correspond to the latter and are in fact 
beats that are generated spontaneously in the hear ventricles. The are in general wider than normal. We 
may determine their width in order to detect them. 

Thus we have created a chart or a database that gives the QRS type for each QRS, referred to by 
its index in the follow-up of the recorded QRS's. 

Normal and supraventricular complex functions are selected using this chart (stage 5-3). A RR 
signal can then be obtained (age 5-4) by calculating the gap between the successive QRS peaks (see figure 
1). 

A transformation into discrete wavelets that is particularly well adapted within the scope of the 
present invention is the transformation into dyadic wavelets. 

In effect, a problem of reversibility is in general posed with regard to wavelet transform. If the 
start of an RR signal is correctly identified, its transformation into discrete wavelets is reversible: we may, 
starting with each transformation, find the initial signal once again through an inversion operation. 

If, on the other hand, the origin of the RR signal is not correctly chosen, the inverse transform 
operation leads to a distortion and does not allow us to recover the initial RR signal. 

Use of the dyadic transform, which is invariant by translation, allows us to solve this problem. 

Information on these dyadic transforms are given in the work by A. Antoniadis and G. 
Oppenheim, Springer Verlag, 1995 and 
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notably in the following articles contained in this work: "The Stationnary [sic] Wavelet Tranform [sic] 
and Some Statistical Applications" G.P. Nason et al., page 281-299 and "Translation Invariant De- 
Noising", R.R. Coifman et al, pages 125-150. 

Figure 6 shows stages for carrying out the detection identification of high-variability zones of the 
RR signal (stage 4-3). 

The coefficients obtained by wavelet analysis are compared to a predetermined threshold, so that 
we retain only the highest-energy coefficients. In this way, this threshold can be carried out on 
breakdown levels j = 2, 4, and 8 (stage 6-1). 

Then, (stage 6-2) we proceed with a filtering of these coefficients, in order to obtain homogeneous 
high-energy zones. 

To do this, we proceed with: 

- a search for the start of each episode or each high-energy zone, 

- an accounting of coefficients in the episode or the zone, 

- a determination of the end of the episode or zone. 

During the search stage for the start of the episode, we count the number coefficients which, 
within a certain window, are greater at the predetermined threshold. When this number is greater than a 
certain predetermined number N s , this is the start of an episode. 

In the accounting phase, we count the number of successive events that are lower than the 
threshold. If this number is greater than a predetermined number, an episode end is detected. 

This stage 6-2 allows a rather precise determination of episodes. 

In some cases, this is nevertheless not sufficient and a greater degree of precision is necessary. 

For this reason, we proceed with a re-setting of the start and the end of episodes (stage 6-3). 

This resetting stage prevents us from missing or allowing single events to pass us by and, 
moreover, allows us to prevent improper identification of the start of an episode in case of a longer 
episode. 

Single events may occur, which are just as important as successive events. For this reason, it may 
be just as important, for the practitioner or cardiologist, to have information available on single events as 
well as on episodes 
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counting several events. Moreover, improper positioning a start of episode, which can, notably, appear as 
a representation on a visualization screen, may trouble this practitioner or doctor and cause him to doubt 
the reliability of the device that is implementing the analysis method according to the invention. 

In order to solve this problem, after identifying the homogeneous high-energy zones of the RR 
signal transformed by wavelets (stage 6-2), we determine the high variability of the numerical signal RR 
(stage 6-3). 

According to one example, we recover the RR numerical signal, and we detect variations in this 
signal, which are greater by a certain predetermined percentage, example 20% (we predetermine a 
threshold of between 10% and 40% or between 15% and 25%). 

When there is such a variation from one point to another in the RR signal, we can say precisely 
that there is a start of an event at this point. 

In the same way, we determine the end of episodes or events. In this way, for each episode, we 
can obtain an interval Nl - N2, which identifies the range of RR signals located between the start and end 
of episodes. 

This episode is then reported on different breakdown levels (figure 7) and only the coefficients 
within this interval are kept as corresponding to the episode. 

Thus it is possible that some coefficients, which are considered after stage 6-2 as being part of the 
start or the end of a homogeneous high-energy zone, are excluded from this zone due to the fact of a slight 
modification, due to resetting, start or end of an episode, or else that coefficients located at the transition 
of homogeneous high-energy and low-energy zones are re-classed from one zone to another. 

Figure 10 shows a patient 20 outfitted with a device 22 of the "Holter" type, with its electrodes 24, 
26, 28, 30, 32 arranged on the patient's body. The Holter records EKG data for a certain period of time, 
for example about 24 hours or over a shorter period (for example a few minutes), or over longer periods 
(for example several days or weeks), and detects and classifies QRS complex functions. 

As illustrated in figure 1 1, the data are, for example, recorded in a memory board 34 of the Flash 
memory type, for example of around 20 Mbytes. Other devices allow us to make recordings on cassettes 
inserted into the Holter 22. 
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Figures 9 and 10 show schematically the essential electronic components of a Holter device. 

Figure 12 shows a device 40 for exploiting the signals recorded on a memory support 34, for 
example a memory board. The support 34 is introduced into a player 36, for example a PCMCIA player, 
and the data are the subjected to handling according to the stages described above. 

Device 40 is therefore programmed to carry out these stages of handling. 

Such a device consists of a micro-computer PC 42 to which the data read by player 36 are 
transmitted through a connection 41. 

More precisely, micro-computer 42 consists of (figure 13) a microprocessor 52, a RAM 54 
memory assembly (for data storage), a ROM 56 memory (for storage of program instructions). 

Optionally, a data acquisition card 60 puts the data supplied by the player into the format required 
by the micro-computer. These various elements are connected by a bus 62. 

Peripheral devices (screen or visualization device 64, mouse 66, keyboard 68) allow an interactive 
dialog with a user. In particular, the means of visualization (screen) 64 allow us to provide a user with a 
visual indication relative to the EKG recorded, to calculated RR intervals, the coefficients calculated by 
wavelet transforms (as in figure 7). By selecting from menus and clicking on buttons or icons, an 
operator, for example a cardiologist, can visualize these different curves, in whole or in part, by enlarging 
them. Set intervals may appear these curves as being episodes of atrial fibrillation, and notably the start 
and end of each episode. 

In micro-computer 42, these data or instructions are loaded in order to implement data handling 
according to the invention. 

These data or instructions for data handling may be transferred into a micro-computer's (42) 
memory zone using a diskette or any other support that can be read by a micro-computer or computer (for 
example, hard disk, ROM dead memory, live dynamic DRAM memory or any other type of RAM 
memory, compact optical disk, magnetic or optic storage element). 

Posting of episodes classified as being atrial fibrillation episodes does not, of itself, constitute a 
diagnosis. Only a cardiologist 
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can, upon seeing the initial electrocardiogram, (which has itself been digitized, but in memory on the 
support 34 and transferred to memory on microcomputer 42), confirm the existence of a pathology. In 
particular, the doctor must still exercise an observation of the base line 3 of the electrocardiogram (see 
figure 1) in order to confirm whether an episode, which was apparently classified as being of the AF type, 
is indeed the sign of an atrial fibrillation. 

The handling process presented above is therefore a diagnostic aid, but is not itself a diagnosis, 
which demands an additional medical interpretation stage. 

EXAMPLE 

This example involves a population of 50 subjects : 19 suffering from chronic atrial fibrillation 
(CAF), 15 from paroxystic fibrillation (PAF), and 16 control patients with a normal sine rhythm. (SR). 

All the subjects were monitored for 24 using a numeric EKG Holter system (DuoHolter, Novacor, 
Paris, France). The recordings were directly digitized at a frequency of 100 Hz on the recorded; they 
were then reviewed and edited manually, if necessary, in order to get a list of consecutive RR intervals. 

Signal handling was carried out with MATLAB® 5.0 software, and statistical analysis of results 
with Statview (4) on a Power Macintosh ®. The results were presented in the mean ± standard-deviation 
form. 

The results were as follows: 

a. EKG recordings of the control group with normal rhythm 

Standard examination of EKG recordings over 24 hours of the 16 SR subjects showed the absence 
of all cardiac arrhythmia. 

In one subject, two AF episodes were falsely detected (duration of 154 and 193 beats respectively; 
mean cardiac frequency of 81.5 ± 14.8 and 76.0 ± 9.5 beats per minute (bpm) respectively). There was no 
significant difference in the temporal indices between the subject with the false detections and the other 
15. 

In this way, specificity in the detection of the AF based on QRS complex functions reaches 
99.96% (1,083,537 Qrs true negatives against 447 false positives). As all the false positives were 
obtained in a single patient, specificity was 
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93.7% if we consider the patient to be a unit of analysis (16 true negatives against 1 false positive). 

b. Group of chronic AFs 

Standard examination of EKG recordings over 24 hours showed permanent AF in the 19 patients. 
The recordings were over a period of 18.8 ± 4.6 hours. In the 19 patients, the number of QRS was 87,953 
± 27,491 and cardiac frequency was 75.4 ±1.0 bpm (interval of 23.to 289.7 beats per minute). 

An SR was falsely identified in 7 patients over 40 episodes (cardiac frequency 82.5 ± 9.5, duration 
374 ± 360 beats); one patient presented 24 of the 40 episodes whereas the others appeared in 6 other 
patients. These false negative episodes were due to atrial flutter and atrial tachycardia in, respectively, 34 
and 1 of 40 cases. In 5 of the 7 patients, the AF was not recognized at the start of the recordings; still, the 
delay in recognition was limited to 178.4 ± 70 beats. 

In this way, detection sensitivity to AF reached 99.2% when we consider the QRS complex 
functions (2,064,197 QRS true positives as opposed to 16018 false negatives). Except for undetected 
episodes, all patients were recognized as suffering from AF and the sensitivity achieved by taking these 
patients into consideration was 100%. The greatest undetected period reached 2,075 beats and the 
shortest 120. The patient with 24 undetected episodes had 10490 beats as false negatives; for this 
particular recording, the detection sensitivity of the AF based on the QRS complex was 89.9%. 

c. Group of intermittent AFs. 

The 15 subjects with PAF had a total 36 PAF episodes, whose length was 2.5 ± 3.4 hours (from 34 
sec to 15h36). During these episodes, the mean cardiac frequency was 86.2 ± 20.7 (from 16.4 to 163.1) 
bpm and the number of QRS was 15,055 ± 22.125 (from 52 to 94,990 QRS). 

Only one of the 36 episodes remained undetected (260 QRS, cardiac frequency from 123.2 ± 28.1 
bpm). For the episodes detected, there were lags or leads in the recognition of the AF at the beginning 
and end of the events: 

- start of episodes: in 8 episodes there was a delay of 91 ± 95 beats (from 2 to 25); 23 episodes 
were detected prematurely (due to supraventricular tachycardiae) of 186 ± 621 beats (from 2 to 
2,994). 

- End of episodes: two episodes were premature (one of 3,486 beats, the other 333); 21 had a 
period of 60 ± 103 beats (from 10 to 429). 
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During the PAF episodes, 14 false interruptions were observed, among which 13 were in the same 
patient, with a duration of 1,484 ±151 beats; the last was 2,001 beats in another patient; the total of QRS 
complex function false negatives then reached 28,695 (1.9%). These false negative episodes were due to 
atrial flutter and to atrial tachycardia in, respectively, 5 and 8 of the 14 cases. 

There were also 37 false detections of 1,573 ± 2,309 beats (from 147 to 1 1,491). 

Likewise in the PAF group, sensitivity based on the QRS complex functions reached 95.8% and 
specificity 91.5%. 

All patients in the group were recognized as having the AF, therefore sensitivity based on the 
patients was 100%. 

Wavelet transform and associated fractal analysis according to the invention furnished a sensitivity 
of 95% for detecting WRS complexes in AF, and a specificity of 99.6%, with a limitation due to delays in 
recognizing the AF at the start and end of episodes of arrhythmia. When based on the patients, specificity 
reached 93% in the normal control group, and sensitivity 0% whether in the CAF or PAF group. 

Intra-atrial recordings were often used to detect the AF, as well as analysis of the electrical activity 
in the wavelet at the level of the esophagus, with a sensitivity and specificity ranging from 52.4-96.6% 
and 92.0-92.3T respectively (see notably S. Rajala et al., "ECG [sic] findings and survival in very old 
people", Eur Heart J. 1985, volume 6, pages 274-252). 

Such recordings were not used given that the goal was the detection of the AF in a non-invasive 
manner with the standard Holter recording. 

For automatic detection, Holter recordings were analyzed using neural networks (T.F. Yang et al., 
"Artificial neural network for the diagnosis of atrial fibrillation", med. & bio. Eng. & Comp. 1994, 32, 
pages 615-619), giving a sensitivity of 92% and a specificity of 92.3% in the recognition of QRS. By 
using a neural network supplied both by information from the RR intervals and a morphological analysis 
of the base line of the EKG tracing, Cubanski D. et al (J. Cardiovasc. Electrophysiol. 1994, 5, pages 602- 
608), working with two consecutive groups of 10 QRS, obtained a sensitivity of 82.4% and a specificity 
of 96.6%. In this way, sensitivity and specificity obtained with the procedure according to the invention 
are better than those of previous studies. 
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Our method has certain limitations. Use of VFC analysis alone is sensitive to everything thing 
that could degrade the VFC, notably supraventricular extrasystolics, supraventricular tachycardiae and 
other types of arrhythmia. Detections of false positives may be due to the presence of a few 
supraventricular arrhythmiae, either isolated or in groups, that generate a significant VFC; detections of 
false positives occurred in 0.04% of the QRS's in the normal control group and in 4.7% of QRS's in the 
FAP group. False QRS's may be correctly reclassified by using a tracing recorded by automatic storage 
of events. Conversely, non-detections may occur in other atrial arrhytmiae, as in paroxystic 
supraventricular tachycardiae, or atrial flutter, due to the relative regularity of these rhythms. 

Under clinical conditions, non-detection is more problematic than detection of false positives, 
because in this case the tracing is not stored for the user, which could prevent any later correction. In the 
cases studied, non-detections were not frequent and represented less than 3% of episodes (1 out of 36) in 
the PAF group and less than 2% of QRS's. Improvements to the algorithm could decrease these non- 
detections, to the detriment, nevertheless, of specificity. 

The invention is therefore adapted to the automatic detection of AF in broad populations using 
standard Holter recordings. This is particularly important in elderly persons, in whom AF, both 
paroxystic and chronic, is frequent and presents a major risk of cerebrovascular incidents. Implemented 
in standard Holter devices or in long-term event recorders, the algorithm or the process according to the 
invention allows us to better bring the AF into evidence in order to take preventive measures before 
serious complications appear. 
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CLAIMS 

1. Process for handling an electrocardiogram RR interval signal, consisting of: 

- detection of high cardiac frequency variability zones in the RR signal, 

- classification of these zones into sine rhythm (SR) and atrial fibrillation (AF) zones. 

2. Process according to claim 1, where detection of high cardiac frequency variability zones is achieved 
by projecting or transforming an RR interval signal into a base or family of discrete wavelets, and 
calculation of coefficients d.(k) giving the projection of the RR signal to each discrete wavelet function. 

3. Process according to claim 2, projection or transformation of the RR interval signal taking place in a 
base or a family of dyadic wavelets. 

4. Process according to one of claims 2 or 3, comprising a thresholding stage of coefficients dj(k) over 
one or more breakdown levels j in order to select the coefficients that are higher than a minimum value. 

5. Process according to claim 4, comprising moreover a step for obtaining high-energy homogeneous 
zones. 

6. Process according to claim 5, where the step for obtaining high-energy homogeneous zones consists of 
filtering threshold coefficients. 

7. Process according to one of claims 4 to 6, comprising, for each high-energy homogeneous zone, a 
search for the start of the zone, an accounting of coefficients in this zone, and a determination of the end 
of the zone. 

8. Process according to claim 7, where the search for the start of the zone comprises a counting of the 
coefficients that are greater than the pre-set threshold. 
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9. Process according to claim 7 or 8, counting of coefficients comprising a countup of the number of 
successive events that are lower than a pre-set threshold. 

10. Process according to claim 9, where an end of episode is detected if the number of events counted is 
greater than a pre-set number. 

11. Process according to one of claims 5 to 10, comprising, moreover, a re-setting or identification stage 
of starts and ends of episodes, with the help of variations in the RR signal. 

12. Process according to claim 1 1, where the coefficients that give the projection of the RR signal are 
retained as being part of a homogeneous high-energy zone only if they are located between a start and an 
end of an episode. 

13. Process according to claim 12, comprising, moreover, a display, on a means of visualization, of at 
least a part of the electrocardiogram and an identification, on this display, of atrial fibrillation zones. 

14. Process according to one of claims 1 to 13, the zones of high cardiac frequency variability are classed 
in SR and AF zones by calculating the fractal exponent of each of these zones and comparison to a 
threshold value H 0 . 

15. Process according to claim 14, where the threshold value Ho is between 0.5 and 0.9. 

16. Process according to claim 15, where the threshold value Ho is between 0.7 and 0.85. 

17. Process according to one of claims 1 to 16, where the electrocardiogram was engaged for a period of 
between a few minutes and several months. 

18. Device for analyzing numerical signals from electrocardiogram RR intervals, consisting of: 
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- a means for classifying these zones as sine rhythm (SR) or atrial fibrillation (AF) zones. 



19. Device according to claim 18, the means for detecting the high cardiac frequency variability, 
comprising a means for projecting or transforming the RR interval signal into a base or a family of 
discrete wavelets, and calculate coefficients d.(k), giving the projection of the RR signal on each discrete 
wavelet function. 

20. Device according to claim 18 or 19, the means for classifying the high variability zones into sine 
rhythm zones and atrial fibrillation zones consisting of a means for calculating the fractal exponent for 
each of these zones and to compare it to a threshold value H 0 . 

21. Device according to one of claims 18 to 20, comprising moreover a means of visualization for 
displaying the atrial fibrillation zones. 

22. Device according to one of claims 18 to 20, comprising moreover a means of visualization for 
visualizing the atrial fibrillation zones and electrocardiograms, and for visualizing the atrial fibrillation 
zones relative to the electrocardiograms. 

23. System for determining cardiac rhythm, consisting of one or multiple recorders for recording 
patients' electrocardiograms, and a device according to one of claims 18 to 22. 
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